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We consider a nearly free falling Earth satellite where atomic wave interferometers are tied to a 
telescope pointing towards a faraway star. They measure the acceleration and the rotation relatively 
to the local inertial frame. 

We calculate the rotation of the telescope due to the aberrations and the deflection of the light in 
;ravitational field of the Earth. We show that the deflection due to the quadrupolar momentum 
e gravity is not negligible if one wants to observe the Lense-Thirring effect of the Earth, 
e consider some perturbation to the ideal device and we discuss the orders of magnitude of the 
e shifts due to the residual tidal gravitational field in the satellite and we exhibit the terms 
which must be taken into account to calculate and interpret the full signal. 

^Within the framework of a geometric model, we calculate the various periodic components of 
tffS\signal which must be analyzed to detect the Lense-Tirring effect. We discuss the results which 
support a reasonable optimism. 

C^ls a conclusion we put forward the necessity of a more complete, realistic and powerful model 
in order to obtain a final conclusion on the theoretical feasibility of the experiment as far as the 
otTServation of the Lense-Thirring effect is involved. 
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I. INTRODUCTION 

The quick development of atomic interferometry during the last ten years 
is impressive. The clocks, the accelerometers and the gyroscopes based on 
this principle are already among the best that have been constructed until 



I i now and further improvements are still expected. This situation favors a 
^jjy renewal in the conception of various experiments, such as the measurement 
of the fine structure constant or the tests of relativistic theories of gravita- 
tion currently develo ped by classical means (gravitational frequency shifts, 
k> , equivalence principle|l6|. Lense-Thirring effect [l7j. etc.. . .). 

The performances of laser cooled atomic devices is limited on Earth by 
the gravity. Therefore further improvements demand that new experiments 
take place in free falling (or nearly free falling) satellites. A laser cooled 
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atomic clock, named PHARAO, will be a part of ACES (Atomic Clock 
Ensemble in Space), an ESA mission on the ISS planned for 2006. Vari- 
ous other experimental possibilities involving " Hyper- precision cold atom 
interferometry in space" are presently considered. They might result in a 
project (called "Hyper") in not too far a future. Most of the modern exper- 
iments display such a high sensitivity that their description must involve 
relativistic gravitation. This is not only true for the experiments which are 
designed to study the gravitation itself but also for any experiment such as 
Hyper where very small perturbations cannot be neglected any longer. 

The present paper is a contribution to the current discussions on the 
feasibility of Hyper. We consider especially the effect of the inertial fields 
and the local gravitational fields in a satellite [l8|. 

There are two kinds of gravitational perturbation. 

1 . The masses in the satellite produce a gravitational field which is not 
negligible. In some experiments, the mass distribution itself can play 
a role : This is, for instance, the case for GPB. However, some other 
experiments are only sensitive to the change of the mass distribution 
with the time. This is the case of Hyper where a signal is recorded 
as a function of the time and analyzed by Fourier methods at a given 
frequency. The modification of the mass distribution is due to me- 
chanical and thermal effects. It depends on the construction of the 
satellite, the damping of the vibrations and the stabilization of the 
temperature. We will not study these effects which can be considered 
as technological perturbations. We do not claim that these perturba- 
tions are easy to cancel but only that it is possible in principle while 
it is impossible for tidal effects from the Earth. 

2. The perturbations due to the gravity of far away bodies (the Earth, 
the Moon, the Sun and the surrounding planets) is the subject of the 
present paper. It is impossible to cancel their action. 

The aim of this paper is to study the gravity in a nearly free falling 
satellite where the tidal effects remain. 

The experimental set-up is tied to a telescope pointing towards a " fixed 
star". However, it experiences a rotation: the so called " Lense-Thirring" 
effect. It has been recently noticed that atomic interferometers display 
a sensitivity high enough to map the gravitomagnetic field of the Earth 
(included the Lense-Thirring effect). This could be one of the goals of the 
Hyper project pj. The effect is so tiny that we will concentrate on this 
question. 

We consider the case where the experimental set up is built out of several 
atomic interferometers similar to those which are currently developed in 
Hannover @ and 0, Paris p| and Orsay [j| and (see section IVI) . 
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In section 1 we introduce the metric, g a p, in the non rotating geocentric 
coordinates and we define a book-keeping of the orders of magnitude. 

In order to study the local gravitational field in the satellite, we chose 
an origin, O, and a tetrad u 1 ?-. which defines the reference frame of the 
observer at point O. The time vector, u? y is the 4- velocity of O. The 
space vector u?U defines the axis of a telescope which points towards a 
"fixed" far away star. The tetrad is spinning around u"^ with the angular 
velocity, vj. 

In section 2, in order to define precisely the tetrad we study the apparent 
direction of the star. 

Then, in section 3, following Q and Q we expand the metric in the 
neighborhood of O (the NiZiLi metric). 

Finally we calculate the response of the experimental set-up and we 
emphasize the interest of spinning the satellite. 

An ASU delivers a phase difference, Sip, between two matter waves. The 
phase difference is the amount of various terms. Some of them can be 
computed with the required accuracy; they produce a phase difference Sipk- 
Then Sip = Sipk + Sip u where Sip is measured. Therefore one can consider 
that the ASU delivers 8ip u . This is this quantity that we want to calculate 
here. 

In this paper, we point out the various contributions to Sip u with their 
order of magnitude. The method that we use to calculate Sip u is a first 
order perturbation method. A more precise method, valid for Sip is now 
available It gives the possibility to model the ASU and therefore to 
study the signal due to the various perturbations which are expected. 

II. GENERALITIES 
A. Conventions and notations 

In the non rotating geocentric frame we introduce the coordinates x a 
with a = 0,1,2,3. We define the time coordinate a; = ct where c is 
the celerity of the light in the vacuum. The space coordinates are x k 
(in this paper the Latin indices run from 1 to 3). We use the notations 
r = {x k } = {x,y,z} and we define the spherical coordinates {r, 9, ip} , 
i.e. x — r sin 9 cos ip , y — rsin#sin</? , z — rcos9. Therefore 
r = || r || = \/ x 2 + y 2 + z 2 . 

We consider an Earth satellite and a point O which is chosen as the 

origin of the local coordinates in the satellite. We assume that the position 

of O is given by its three space coordinates, r = {x, y, z} — {x k } , as three 

known functions of the coordinate time, t. Then we define the velocity of 

_ _> dr 
point U as v = — . 

dt 
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The proper time at point O is s. The motion of O can be described 
as well by the four functions x a = x a (s) . The four-velocity is defined as 

a 

U = ——. 

as 

The Newtonian constant of gravitation is G. 

We use geometrical units where the numerical value of G and c is equal 
to 1. 

The metric tensor is g a p\ it inverse is g Q/3 . The Minkowski tensor is 
T) a p = diag [1, — 1, — 1, — 1]; its inverse is i] al3 . 

We use the summation rule on repeated indices (one up and one down). 

d( ) 

The partial derivative of ( ) will be noted ( ) = d a ( ) = — — . 

ox a 

In the sequel we introduce different tetrads, i.e. a set of four vectors ej? 
labeled by the means of an index called a or (a) , such as e%g a p,e^ = r\%p. 

The "Einstein" indices, a, j3, a, etc., and the "Minkowski" indices a, (3, 
etc., can be changed one into the other by the means of the tetrad : ( ). = 
e & ( )u • ^he me t r i c tensor is used to raise (or lower) the Einstein indices 
while the Minkowski tensor is used for Minkowski indices. 



B. The fundamental element 



In the sequel we consider the Parametrized Post Newtonian theories 
jlCj . The relevant PPN parameters which appear below are 7 and a\. The 
parameter 7 is the usual parameter connected to the deflection of a light 
ray by a central mass. The parameter a\ couples the metric to the speed, 
—w, of the preferred frame (if any) relatively to the geocentric frame. In 
general relativity, a-y = and 7=1. 

The preferred frame is generally considered to be the rest frame of the 
Universe where the background radiation is isotropic (||to|| ~ 10~ 3 in geo- 
metrical units). 

Let us define now several quantities which will be used in the sequel : 

• 2Mq is the Schwarzschild's radius of the Earth (M ffi ~ 4.4 mm). 
As we use geometrical units (G = c = 1), M@ is also called the "mass" of 
the Earth. 

• J® is the angular momentum of the Earth in geometrical units. 

1 + 7 + Oil /4 - 

The relevant quantity which appears below, is J = Jq. We 



define J = 



= J m ~ 145 cm 2 . 



J A r 1 

g = —2 — ^ 1 — a\U w is the definition of g. 

r A 2 
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• U is the Newtonian potential 

C / = ^(i-j 2 (^) 2 p 2 + a) +C /, (1) 

where is the radius of the Earth and £/* the potential due to the Moon, 
the Sun and the planets^^- I n spherical coordinates the Legendre poly- 
nomial P 2 reads P 2 = - (3 cos 2 9 — l) . The quadrupole coefficient J 2 is of 

order of 10~ 3 and A represents the higher harmonics; it is of order of 10 -6 
[Tl|. It depends on the angle tp and on the time t because of the rotation 
of the Earth. 

In the non rotating geocentric coordinates the significant fundamental 
element is 

ds 2 = (1 - 2U) dt 2 + 2g ok dx k dt - (1 + 2jU) S jk dx } dx k (2) 

where (g) k = — (g) k = gok- In the expression (J2J, we have dropped post 
Newtonian corrections which are too small to be considered here. 



C. Orders of magnitude 



We consider a nearly free falling, Earth satellite on a nearly circular orbit 
of radius r (with r ~ 7000 km) . 



1. Orbital data. 



The velocity of the satellite is of order 0\ ~ J (i.e. 0\ ~ 2.5 10 5 
for r ~ 7000 km). 

In the neighborhood of the satellite the potential U is of the order of 

the potential of the Earth U ~ — — ~ 2 ; the coefficients g fc fulfills the 

2 r 

relation \g ok \ > ^ f^j ~ 750 4 » 4 (M)- 

The post Newtonian terms that we have neglected in the metric (JSJ are 
of order O4. 



2. The size of the set-up. 

Now we define d ~ rO\ (i.e. d ~ 175m for r ~ 7000km). We assume 
that the size of the experimental set-up in the satellite is X ~ e <i where 
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e < 1. In the sequel we will consider an atomic Sagnac unit the size of 
which does not exceed 0.9m (then e ~ 5 10~ 3 ). 



3. The local acceleration. 



The acceleration, || a || , can be measured by accelerometers comoving 
with the satellite. It is called "the acceleration of the satellite relatively 
to a local inertial frame". The satellite is nearly free falling, therefore 

I|~a1| << a < 2 x -. 

r 

Our assumptions are summarized below. The numerical values are ob- 
tained for the Earth with r ~ 7000 km, M m /r ~ 2 ~ 610~ 10 and 
v/c ~ 0\ ~ 2.5 10~ 5 (i.e. u ~ 7.5km- s^ 1 ) where v is the velocity of 
the satellite relatively to the geocentric frame: 



M e /r - 2 ~ 



X < eOir = e<i, with e < 1 



|| ~a\\ < a ~ 2 x 



M® ~ 4.4 mm 


J ~ J® ~ 145 cm 2 


r ~ 7 000 km 


-4- - 750 » 1 



Oi ~ d/r 


O2 ~ fa 


| ff ofc| - 750 O4 


4 


2.510 _b 


610- i0 


310~ ib 


410- ly 



<i ~ 175 m I a ~ 8 m • s 2 



Table 1 

In the sequel we will assume the preceding relations. 



D. The NiZiLi comoving metric 



Before closing this section dedicated to generalities, we give the expres- 
sion of the fundamental element associated to the NiZiLi metric. 

First we choose an origin, O, in the satellite and we choose a tetrad, e% 
whose vector e? = u a is the 4- velocity, u a , of point O. Thus, the vectors 

define the basis of the space vectors for the observer O. The coordinates 

which are associated to the tetrad are the space coordinates X k and the 
time X" = T. 

Following the procedure defined in Q and @ we find the fundamental 
element and the local metric tensor 
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ds 2 = G & 0dX & dXP +e 2 6 with (3) 

- R m xkx " 3 \ R m.i x ' kx ' 3xi + ••■ 

2 1 - 

G 0m = ^rhk X - z R 0krh3 X X ° ~ 1 R bkrh' 3 ,l X X3 X + - 
Gnm = Vfim- 2 R nkAj XkXJ ~ Q R nkm3A XkXJX£ + ■■■ 

where we have used vector notations i.e. V for j V k j , V ■ W for V * 14 7 * , 

etc. Every quantity, except the space coordinates X k , are calculated at 
point O. Thus they are functions of the time T. 

R &i3ap, 1S tne R-iemann tensor obtained from R a a a ^ at point O : 

RaBafi = r« — Bfi.a ^a — 3a.fi ~\~ ^ /3a^s — afj, ^ (3/jJ^s—cta (4) 

where r a _0 M is the Christoffcl symbol. 
Q~~ k is the antisymmetric quantity 

de" de 13 \ \ 
Vf k = £ ( %,.k - 96k, ) +^ll e iX-X e i W ( 5 ) 



"j* 2 V y °i' fc aofe ^7o 2 U 3 ds ds * 

\\ / / o 

Due to the antisymmetry of the quantity X k dX m which 

is present in the expression of ds 2 can be written as fK^ X k dX m = 
• dX. The space vector fl is the physical angular velocity. 
It is measured by gyroscopes tied to the three space orthonormal vectors 

The vector a is the physical acceleration which can be measured by an 
accelerometer comoving with O. It is the spatial projection at point O of 
the 4-acceleration of point O. 

At point O (i.e. X = ) the time T is the proper time delivered by an 
ideal clock comoving with O. 

The first tetrad that we consider is actually called e% : 
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e? = u° = 1 + — + U + 4 , eg = u fe = \1 + — + Uj v k + 4 
e\ = U+^ + U^v k + 1 Uv k -g 0k + O± (6) 

Calculating Hq one finds 



o = n LT + n dS + n Th (7) 



k 




n LT ) ~ l---(j.ft) ft -^-VUAw] (8) 



(fids)* ~ ((l+ 7 )VC/A t /) fc and (o T ,)^ Q«A |j (9) 

Qlt is the Lense-Thirring angular velocity, Qds and $7^^ are the de Sitter 
and the Thomas terms [2l| . 

Now it is straightforward to calculate the order of magnitude of the 

J c 

Lense-Thirring angular velocity. One finds — 2O4- ~ lCP 14 rad • s _1 . 

We will see that it is relevant to limit the expansion of the metric at 
order e 2 Oq; therefore we consider only the linear expression of the Riemann 

tensor above and we neglect he term ( a ■ X ) in the metric [21 



fc 



In the sequel we consider an other tetrad u"^ . Except for the notations, 
the general results above still hold with the associated coordinates. 



III. ABERRATION AND DEFLECTION OF THE LIGHT 

In the satellite, the experimental set-up is tied to a telescope which 
points towards a "fixed" star. We assume that the star is far enough for 
the parallax to be negligible. However it is necessary to account for the 
gravitational deflection of the light ray and for the aberrations in order to 
describe the rotational motion of the telescope during the orbital motion 
of the satellite. 

In space time, the direction of the light from the star is given by the 
f 9^ cp 1 

4- vector L a = < 1, — — > where ip is the phase of the light. In order to 
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calculate the phase tp(t,x k ) at point {x k } and time t we use the method 
which is summarized in annex A |12| . 

The monopolar term of the Newtonian potential gives 

where £ is the unitary vector of figure ^and ~n —~r /r . 
Several contributions must be added to the expression ifTDl) . 

L a = Sl, -l k + (1 + 7) M.HL^L + 5i k + SL k 
I r i-n.£ 

• The term S£ is due to the quadrupolar term of the Earth. We give 

its value in appendix A. Its order of magnitude is JiOi ~ 10 _3 C>2 when £ 
is nearly orthogonal to the plane of the orbit. 

• The term due to —ctiUw k , a part of gok in the metric @, results 

(cx\ w • £ \ 
1 r . The demonstration 
2(1 + 7); 

is straight forward. 

• The contribution due to the rotation of the Earth has already been 
considered in the literature f[l^[). Since it is proportional to the ratio of 
the angular momentum and the square of r, it is negligible. The same 
conclusion holds for the term A in expression Q. 

• The Sun, the Moon and the other planets, give a contribution due 
to {/» in QJ; it varies slowly with the time and it is negligible, especially 
within the framework of a Fourier analysis at a much higher frequency. 
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For the observer O, the space direction of the light is the four vector 
A Q = L a — Lp vP u a . The components of X a relatively to the tetrad e& are 



X& = e^X/i 



-h 



-L a e% 



-Aa 







In order to calculate X k , we use the expression (JBJ for the tetrad. We 
normalize A, i.e. we define A = A A such as — A a A Q = A • A = 1. The 
tetrad © is especially useful to catch the orders of magnitude of the various 
terms involved. However it is not the comoving tetrad that we are looking 
for because the telescope that points towards the far away star rotates 
relatively to this tetrad. The angular velocity of the telescope relatively to 



HI is 



f2* = A A — . Straightforward calculations give 
at 



A fc = 



i+e-v- -v 2 



+ (1 + 7) 



M n-l 
r \-n-i 



i- ( si- 



■5L) 



lk M 1 + 7 
r 1-H-e 



6£ + 5L + 3 



and 




+ (tAvf 



t-v - 



(11) 



■I 



Let us notice that we neglect the terms of order — x O4, which are much 

r 

smaller than the Lense-Thirring angular velocity f~ 750 O4 x — ). 

r 

Now we introduce the tetrad uf , which is obtained from through a 



pure space rotation (i.e. it 



(0) 



and whose vector points 
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towards the far away star (it"^ = —A"). The rotation of the tetrad |tt^ j 
relatively to | e ^| is characterized by the most general angular velocity 

f2„/ e = 17* — tuA where —zuA is an arbitrary angular velocity around the 
apparent direction of the star. 

IV. THE LOCAL METRIC IN THE SATELLITE 
A. The relevant terms 



Now we need to calculate the NiZiLi metric associated to the tetrad u? <, . 
We limit the accuracy of our development to a few tens per cent in order 
to take advantage of important simplifications. 

Let us assume that the fundamental element is known : 

ds 2 - (1 + K (om ) dT 2 + 2K mk) dTdX^ + ( V{k)U) + K (k)U) ) dX^dX® 

We want to calculate the effect of Kt a \{p-\ on the atom interferometer in 
the satellite. In order to use the method summarized in the appendix A we 
calculate the quantity W which gives the physical quantities that one can 
measure : 

* = K mo) + 2K mh) vj?> + K {k){j) v^v^ (12) 
where Vg is the velocity of the atoms (i.e. the unperturbed group veloc- 

ity). 

We are dealing with slow cold atoms, therefore v g ~ r\0\ with r\ << 1 
(we can take the value rj ~ 0\ valid for current experiments). In the 

Lense-Thirring term arises from 2K^^v g k ^ ; it is of order of 2 — Xv g ~ 

1500 r\eO§. If the expansion of K(o)(o) is limited to order £ 2 Oq, the accuracy 

is of order ~ 13% with v a ~ 20 cm • s . 

1500r; 9 

In the case that we consider, the following orders of magnitude hold true 
: £ ~ \/Oi — 0\li an d f] ~ 0\. Therefore, it is not necessary to consider 
terms smaller than O7 ~ e 2 Oq in J^(o)(o) 7 smaller than O5 in Kr \ijA and 
smaller than O3 in K(u\ff, . The same relations hold in the coordinate system 

associated to the tetrad {"fcr)} as weu as { e !s} b ecause both are next to 
each other. 

It can be proved that in the expression ^ the time can be considered as a 
parameter and that, within the present approximations, it is not necessary 
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to consider time derivative in the expression [3J Therefore one obtains 
Ggg = 1 + 2 a ■ X - U, t3 X k X$ - ^U, m X % X*X l + e 2 6 

( -. -» -\ rh 

Gg A = — |fio A Xj + e 2 5 and G ftlh = r] flfn + e 2 4 

where S7o is given above (see expression 0) while the expressions such as 
U are nothing but (u^mn e ]f ■ The position of the observer changes 
with time, therefore this quantity is a function of T. 

In Ggg, such an expansion is limited to the terms of order of e 2 Oe, thus 
the approximation e™ — 5™ is valid and U,j,.^ (U.kj)o • For the same 

holds true : (JIq^J ~ (tio^J ( se e expression^ . Therefore, one can identify 

the space vectors ~e ^ of the tetrad and the space vectors d k of the natural 
basis associated to the geocentric coordinates. This would not be valid 
with an higher accuracy where terms smaller than e 2 0§ are considered. 

The change of the tetrad u? a -. < — > ej? is just an ordinary change of basis 
in the space of the observer O. In this transformation, dT, Ggg = G( )(o), 
G 6lh dX A = G mk) dX^ and G nA dX A dX n = G (])(k) dX^X^ behave 
as scalars. 

As a consequence, we put forward that within the required accuracy, 
it is possible to give a simple description of the Hyper project with the 
Newtonian concept of space. It is now straightforward to calculate ^ : 

vf = 2a-X-U, (k)u) XWXW -±U, (k)um X<MX<»XW 

{ (fio + d*) A x} (k) + e 2 O e 

(fe) 

Now we assume that any quantity can be known with an accuracy fi of 
order of 10~ 5 < 0\. It means that the geometry of the experimental device 
is known with an accuracy of order of 10/Ltm. The position of the point O 
on the trajectory of the satellite is known with an accuracy better than 
70 m, the eccentricity can be controlled to be smaller than 0\, etc. 

We consider that 'J is the amount of two terms, and ^f u : the term 
can be modelled with the required accuracy while ^ u is unknown. The 
terms VE^ fulfills the condition /ix^k < e 2 Oq. With the previous orders of 
magnitude (e ~ 0\ii ~ 5 10~ 3 and [i^ 0\ ~ 2.5 10 -5 ) one finds 



* u = 2a-X-U, (k)u) I%( J »-2{l]Ai}^ 9 +£ 2 6 /^ 
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FIG. 2: The satellite and the fixed star. 



U. 



0)(i) 



is obtained from the expression of U limited to the monopolar 



and the quadrupolar terms and 



Ult-ihA-AA— + — 
at I r 



where Qlt is obtained from [S] and 



d v 
~~dt 



a - VU: 



B. Calculation of \P„ 

We consider that the motion of the satellite is aNewtonian motion which 
takes place in the (x, y)-plane while the vector I lies in the (x, z)-plane. 
We assume that the eccentricity, e, does not exceed 0\. 

We define 



J = J x e 1 + J y e. 2 + J z e§ , n — cos ej + sin eg 
—I = cos a e j + sin a eg , w = w x ej + w y eg + iu z eg 

J, £ and wJ are constant vectors. The angle and the distance r depend on 
the time T. 

First we define the special triad u'^ : 



= cosa e j + sma e g, = sma e j - cosa e §, u, 3 n 



14 



Let us outline that we have defined (t? (1) ) = - (a) + Oi. Then, 
in order to obtain the final tetrad m"^, we perform an arbitrary rotation 
around A : 



«(i) 
w (3 ) 



-A : 



"(1) 



u', 2 s cos a + u',^s sin cr + 0\ 
—u'^ sin cr + u'^ cos ct + 0\ 



(13) 



(icr 



where — — = — w is the angular velocity of the triad relatively to 

{*«}■ 

We can now assume that the experimental set-up is comoving with the 
triad ui n \ whose vector ui\\ points towards the fixed star. 

In order to calculate the quadrupolar term in U, we define the unitary 

vector, k along the axi-symmetry axis. One obtains 



M ffi 1 M e 

u — -Ji — 

r 2 r 



Therefore reads 



with 

Q 




(Lit) 



a\ Mj 
4r 2 



n A w Ax 



—2m (m(i) Ax) • w ff 
-2—^ ((«(i) A n) A x) ■ v g + 2 A a) A f) • v g 

+— ( x - 3 (n • .t) J + - J 2 — ( — J Q 
+2a-x 



.4 

I? 

C 

D 



(14) 



1 — 5 ("fc ■ t?) j A 2 - 20 (~fc • 7?) (r? • a) ( ~k ■ x) 
+5 (j (y -7?) 2 -3^ ^-x) 2 + 2("fc -x) 2 



and 
~J = J~k 



In \I/ U , the two terms 2a ■ X and — 2tn ^U(i) A • w g are directly cal- 
culated in the local coordinates associated to the tetrad u"^ moreover, to 
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zone (1) (2) (3) 
Lasers JJ, JJ, j| X«n JJ, jj, j| 



^ ft tt tt **"tf ft tt 

ay$^it D *»* i i ( 3 > 

FIG. 3: An atomic Sagnac unit (ASU). 



calculate the other quantities where n is involved, one can drop the terms 
of order 0\ in the expression of W(i). The reason is that the quantities that 
we drop are either negligible or included in ~^>u- Therefore, it is clear that 
we can consider the space as the ordinary space of Newtonian physics and 
that the usual formulae to change the basis d k into e j or ~u (/^are valid. 

In i|14fl . the terms of lines A, B and C are due to various rotations : 
respectively the Lense-Thirring rotation, the spin around the view line of 
the star and the aberration. The term of line D gives the gravitational 
tidal effects (which are mainly due to the Earth) and the last term (line 
E) corresponds to some residual acceleration due to the fact that point O 
is not exactly in free fall. 

It is now possible to calculate explicitly ^ u with the coordinates X^ k \ 
comoving with the experimental set-up [22] • 



V. THE EXPERIMENTAL SET-UP 



A. The atomic Sagnac unit 

An atomic Sagnac unit (ASU) is made of two counter-propagating atom 
interferometers which discriminate between rotation and acceleration (see 
figure EJa). 

Each interferometer is a so-called Ramsey-Borde interferometer with a 
Mach-Zehnder geometry ( figure yJ-6) . The atomic beam from a magneto- 
optical trap interacts three times with a laser field. In the first interaction 
zone the atomic beam is split coherently, by a Raman effect, into two beams 
which are redirected and recombined in the second and the third interaction 
zone. 

The mass of the atom depends on its internal state, therefore it is not a 
constant along the different paths. However, the change of the mass is very 
small; it leads to negligible corrections on the main effects which is already 
very small. Considering cesium, we assume that the mass of the atom is a 
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constant m = 133 x 1.66 x 1CT 27 = 2.2 x 1CT 25 kg. 

In this case the wave length of the lasers is A = 850 nm. The momentum 

Airh 

transferred to the atom during the interaction is — — . The recoil of the atom 

A 

results in a Sagnac loop which permits to measure the angular velocity of 
the set-up relatively to a local inertial frame. The device is also sensitive 
to the accelerations. 

One can easily imagine that there are many difficulties to overcome if 
the Lense-Thirring effect is to be observed. In particular, the geometrical 
constraints appear to be crucial. In an ideal set-up the two interferometers 
are identical coplanar parallelograms with their center Os and 0' s at the 
same point. One can consider several perturbation to this geometrical 
scheme. 

1. A shift : Os and 0' s are no longer at the same point; 

2. A tilt : The plane of the two interferometers are now different; 

3. A deformations : the interferometers are no longer identical lozenges, 
they are no longer parallelogram, even not plane interferometers. 

The geometry of the device is fully determine by the interaction be- 
tween the initial atomic beam and the lasers ; moreover the geometrical 
description is already an idealized model ; Therefore a full treatment of 
the atom-laser interaction in a gravitational field is obviously necessary to 
study the response of the ASU (see 0). However the geometrical model 
is useful to give a physical intuition of the phenomena. In this context we 
study here the effect of a shift on the signal. 

In the sequel we consider an ideal set-up where the beams are perfectly 
coherent and perfectly parallel, point Os being a perfect center of symmetry 
for the atomic paths. We assume also that the two counter-propagating 
atom interferometers are identical and located in the same plane but that 
the separation of their center of symmetry, Os and Sl is the vector 5 = 

OsO?. 

We will assume that the velocity, v g , of the atoms is of the order of 
20 cm • s _1 and that the size of the ASU is of order of 60 cm, therefore the 
drift time is 2Tb ~ 3 s. 

B. The phase differences 

The configuration which is presently considered in the Hyper project is 
on figure 0] 

We consider one among the four interferometers of the device above. 
Using the results of appendix A, we obtain the phase difference that we 



17 



_^ t> f* telescope 




(!) 



FIG. 4: The experimental set-up. 



want to measure: 



5<p=~ [ ^(t')dt'-^f ^(t')dt' (15) 

z Jt-2T D /{2) 1 Jt-2To/(l) 



The integrals are performed along path (2) and (1) of figure 0-6. The 

2 

Tfi C 

" angular frequency" u> is defined as — - — . 

h 

We consider that the origin is the center of symmetry O5. The coordi- 
nates are X = {X^} . We define v g = —j- as the unperturbed group 
velocity of the atoms 

We make use of the properties of the Newtonian trajectory : 

r = — -, ft = cosOe; + sinflek, = fisr (1 + ecos6*) 2 

1 + e cos dT 

with ils = — tt. We consider the case where the ellipticity, e, is much 
smaller than unity. Therefore, up to first order relatively to e, one obtains 

- = — (1 + e cos 6) , 6 ~ n s t + O + e sin (Q s t + 6 ) 
r r 

We assume e < 0±, therefore in the expression (|14|) of ^ u one can assume 
that r — ro is a constant because the corrections due to the eccentricity 
are included in ^k- 

During the flight time of an atom, one cannot consider that it remains 
a constant in expression l|14l) . One can consider that the coordinate, X = 
X^- 1 ' of the atom is a function of the time: X = v g (T — To) . Then it = 
no+6 n where Tiq depends on T only and 5n on T and X/v g . Therefore, 
in expression (|15fl . we expand the different terms relatively to X/v g . Any 
of the correction can be included in except the first order correction, 
Srpjj, of the first term in line D. 

6M® /_> lA _ _> 

oipD = 3— [ x ■ — ( x ■ n J x 

r 4 V v n 
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The result of the integration is now straightforward. We change the nota- 
tion : Tlq — > n and T — ► T. In usual units one finds : 




where S — v q Tjj is the area of the Sagnac loop. 

Am 

Line B in i|14fl gives no contribution because X , uri) and v g are in the 
same plane. The quadratic terms have disappeared because Os is a center 
of symmetry. 

The two interferometers of the same ASU are assumed to lie in the same 
plane but not necessarily with their center of symmetry Os and 0' s at 
the same point. Therefore adding and subtracting the phase differences 
delivered by the two interferometers one finds the two basic quantities 
which are measured by the set-up i.e. : 



Mi = ^ \<><P -o<p) = -£ [cT D ) <j — it(a) • n + w (3) ■ a + il LT • u (2) 

\2 



2 v r r ' A 
2ir(cT D ) 2 f (a 0s a Q 



c 



,2 



■«(3) 



mc (cT D fs M s 
~2h r 5 r~ 



• ft) (u {3 ) ■ ^) + (t2(3) ■ n) ■ 



1 fx ' , x \ 4tt(cTd) 2 J "a ^ 
M2 = ^ W + M = x \ "^2 ' "(3) 

where now ~a — ao s + ■ In the expression above, we have dropped 

, mc c (cTd) 2 M e 47rft . . 
the term —ttt-j 5 t \ u (3) ' n J ( u (3) ' v ) because it can be 

included in ^>k- 

a Os u< (k) IS 4-acceleration of point Os- Therefore, one calculates 

(ao s — clo') ' "(3) — ~^,(k)(S) x Os wnere we have dropped the non sig- 
nificant terms which can be included in 
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The quantities which can be measured are 

Hi + 2vgfi 2 = — (cT D ) < — 2"U( 3 ) • n + S2 L t • "(2) 

27r(cT D ) 2 M e / t _ ^ . r T\\ 
^ — (u (3) - — -3(n-« (3) ) I //■— ! I 

-f (cT,) 2 ^ 



r V r 

2 



| (t? ( i) • n) (u (3) • + (5(3) • n) • ^ 



1 ^ ' -L x \ ^(cT D f f ^ 
where <5 = ~x ' ~ ~%o s ■ 



Discussion 



We define the projection, J „ of J on the plane of the orbit : J „ = 
J M (cos 9 j ~e i + sin Oj ~e ' 2 ) and to M = io M (cos^ ~e ^ + sin9 w ~e ^) ■ Going 
back to usual units (with M e ~ 4.4 mm and J ~ 145 cm 2 ) we obtain 



/il + 2 — ^2 

c 



2tt (cT d ) 
Xr 



{K + K a + K 2a + K 20 



+Kg_ a + Kg +a + K 2 0- a + K 2 g +a + K 2S - 2a + K 2 g +2a ] 



with 



K 2a 
K 20 



Me ,„ • 2 ,n <5 (3) 

- — 3 sin a — 1 x 

4r v ; r 



v I J 

— ' [(1 — sin a) cos (a + 0,/) - (1 + sin a) cos (<r — 6j)\ x — ^ 



-4 cos a cos a x — 



3M e <$« 

— — cos a sin a sin a x 

2r r 



3M® 
4r 
3M 



(l — sin 2 a) 



^(2) ^(3) 

sin (2a) x h cos (2a) x 



5(3) 

(l-sin 2 a) cos (29) x 

4r v 1 r 
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K0 +a 



M®^ x {-2(1 + sin a) sin (6 - a) 

ro c 

+ — cos a sin (0 — a — 6 W ) x — 
2 c 

/- . \ . //, \ w 3 I 
+ — (1 + since) sm(0 - er) x — > 
2 c J 

x {-2(1 -sin a) sin (6 + a) 

r c 

a\ . . . tu M 

+ — cos asin(f + (T- 0^) x — 

A C 



Ctl ,, . . .10 

— — (1 — sina) sin (0 + a) x — 

^ C 



. U, 



3-2- (1 + sin a) cos (20 - a - 6j) x 



c 

3M 



6 W 

— — cos a ( 1 + sin a) sin (26 — a) x 

4r r 



/Me 



3/2 



CV 9 T D 

2ri 



(1 + sina) cos (26 — er) 



w J 
# 2(?+CT = 3-^ (1- sin a) cos (20 + a - 6j) x 

c r 

+ cos a (1 - sina) sin (26 + a) x - — 

4r r 

M ffi \ 3/2 cw„T n 2 , , n 

9 (l-sina)cos(20 + ( r) 



2r 2 



3M ffi o f 5 (2) (5( 3 M 

^20-2<r = -r — — (1 + sin a) <^ sin (26 - 2a) x cos (20 - 2<r) x \ 

8r { r r J 

i^2e+ 2( T = -^^(1- sina) 2 (sin (2<9 + 2cr) x — + cos (20 + 2cr) x — 1 
8r I r r J 

Each of these terms, except K , has a specific frequency. These terms 
can be measured and distinguished from each other. 

The Lense-Thirring effect due to the angular momentum of the Earth ap- 
pears in the terms K a and K29± a while the possible existence of a preferred 
frame appears in Kg± a which depends on the components of a{w . 

The signal due to the Lense-Thirring effect is associated with the signal 
due to The two signals display the same order of magnitude when 
~ 1 nm. Today, it seems impossible to achieve such a precision, this is 
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the reason why 5^ should be calculated from the Fourier analysis of the 
signal itself, altogether with the angular momentum of the Earth, J and 
the velocity a{w . 

If the sensitivity to measure the Lense-Thirring effect with an accuracy 
of 20% is achieved, it should be possible to know a{w with an accuracy 
better than 1CP 7 . Considering that ~w is the velocity of the rest frame of 
the Universe (\\vj || ~ 10~ 3 ) it would give precision onaj of order of 10~ 4 . 

The interest of the spin is obvious. If a = cte (no spin) the signal is 
the sum of two periodic signals with frequency vo and 2vo where vo is the 
orbital frequency of the satellite ; therefore one ASU gives two informations 
(two functions of the time). When the satellite spins, we get 9 functions of 
the time t. The information is much more important in this case. 

VI. CONCLUSION 

In this paper we have sketched a method to take into account the residual 
gravitation in a nearly free falling satellite, namely the tidal and higher 
order effects. We have shown that these effects arc not negligible in highly 
accurate experiments. 

We have shown that many perturbations must be considered if one wants 
to observe the Lense-Thirring effect and we have exhibited the various 
terms that one needs to calculate in order to obtain the full signal. 

Compared with GPB, the principle of the measure is not the same, the 
difficulties are quite different but the job is not easier. For instance, con- 
sidering the quantities K a or K2e±<j above, one can check that 5^ must 
remain smaller than 2 nm for the corresponding signal to remain smaller 
than the Lense-Thirring one. It does not seem that such a precision can 
be controlled in the construction of the experimental device itself. It is 
therefore necessary to measure d^ 1 ' with such an accuracy. 

In the problem that we have considered, there are 9 unknown parame- 
ters : i) the three component of J , ii) the three components of a{w and 

iii) the three components of <5 . On the other hand, the experimental set-up 
displays 9 periodic functions but the distribution of the unknown param- 
eters among the 9 functions happen in such way that the four parameters 
J and are present in the 3 functions K a or K 2 g± a and the three pa- 
rameters a\w in the two functions Kg± a . Only the two parameters 5^ 
and 5^ are over determined by the four functions K 2 g, K 2lJ and K 2 $± 2a . 

Let us assume that 9 and o are known function of the time (frequency and 
phase). This implies that in the geometric scheme that we have explored, 
one can determine 18 unknown parameters. Therefore J , S and a{w can 
be known and the Lense-Thirring effect can be observed with an accuracy 
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of a few tens percent. The same sensitivity on the phase difference of 
matter waves in the interferometers yields an accuracy of 10~ 7 on a\W 
which would increase our knowledge on a,\ by one order of magnitude. This 
optimistic conclusion must be tempered with the remark that only the shift 
has been considered here while several other geometrical perturbations play 
their role. Moreover a crucial point is the knowledge of the phase of the 
various periodic functions K. The geometric scheme fails to describe the 
change of the phase of the atomic wave when it goes through the laser 
beam and we believe that the preceding conclusion holds only in the case 
where the change of the phase along the two paths differs by a constant. 

As a conclusion, we put forward that only a more powerful model can 
answer the question of the theoretical feasibility. This model should take 
into account all the gravitational perturbations that we have outlined here 
and it should consider the interaction between laser fields and matter waves 
in more a realistic manner. 



APPENDIX A: THE GRAVITATIONAL PHASE SHIFTS 

Let us assume that space time is quasi Minkowskian. Therefore, the 
metric is: ds 2 = {r] a p + h a p) dx a dx$ with \h a p\ << 1. 

Let us consider, at the eikonal approximation, a wave which propagates 
from the point A to a point B. The phase at point A is known. It is 
y>A(t) = u>t where u> is a constant and t = x°/c the time. The phase, 
ipB(t), at point B is the amount of the unperturbed phase, tp B (t), and the 
perturbation dips due to the term h a p in the metric. 

Two different cases are relevant for the problem that we study 

1. Point A is a far away, fixed, star. It is the source of a light wave. 
From the knowledge of at every point B one can deduce the 
apparent direction of the star. Then it becomes possible to determine 
the vector —L a which points towards A. 

2. Point A is the source of an atomic wave which enters a matter wave 
interferometer. At point B interferences are observed on a detector 
(figure EO-fr). The phases <pBi(t) and <pB2(t), of the waves which in- 
terfere at point B depend on the path, (1) or (2), that each wave has 
followed. The response of the detector at point B depends on the 
difference <pBi(t) — fB2(t). Therefore one can obtain the response of 
the interferometer once the phases <fBi(t) and (fB2(t) are known. 

In order to calculate the phase <fiB(t) at point B and time t we use the 
method developed in [l^. We summarize briefly the method for particles 
of mass m (for the light m = 0). 
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First we neglect the perturbation and we consider a point M which moves 
at the group velocity, v k , and arrives at point B at time t. The worldlinc 
of M is x k = x k (t') with v k = dx k /dt' and x k (t) = x%. The point M has 
left A at time tA such as x k (tA) = x k A . The time tA is a function of t. 

Now we dchnc * = h o + 2hokVg+hkjV k v J g where \1/ is calculated at point 
M. Therefore ^ is a function of t' . 

One can prove that the perturbation dips is 

5<p B = % [\(t')dt> 

1 Jt A 

where h u is the energy of the particle. 



APPENDIX B: DEFLECTION OF THE LIGHT DUE TO THE 
QUADRUPOLAR TERMS OF THE EARTH 

The Newtonian potential to be considered is 

where P 2 = ^ (3 cos 2 6 — 1). The unitary vector k defines the axi-symctry 
axis and cos 8 = ^L_L. The quadrupole contribution is due to the term 



r \ r 




Sip = J 2 R t 



where XM{t') = ~r' = xb + s£. Here s is t' — t, it varies from sa = —00 
(when A is far away) to s B = at point £>. 

We calculate the gradient of dip at point B the coordinates of which are 
V = {x k } : 

VJp = -~^J 2 RqQ with 



J —00 1 J — 00 1 J — 00 ' 
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Each of these integral can be exactly computed. The deviation of the light 
ray is A£ — — \75<p. We can eliminate some components parallel to /which 
would in any case disappear in the normalization procedure. One finds 



S£ = 



M e (Rq\ ak + bn 
— — J 2 — I —with 

4r V r 



2 



= 2 



b = 



(it • (it - l>y 
-6 ( it ■ ~k ) ( n ■ 7) + 2 (j ■ ~fc ) (it ■ 
+4 ( n 1 • - 2 (7 • ~fc ) 
-3 (r? • "fc) 2 (7? • T) 2 + 9 (it ■ (it ■ T) 
-2 (r? • ~k) (j ■ ~k ) ("n • 7) + ("n • lQ 2 - 8 • ~fc 
-2 (7 • T) 2 + 6 (it ■ ~k) (J ■ ~k ) - 3 ("n • 7) + 2 



where r = Vr^ and n = rjr. 

When i is orthogonal to the orbit, the maximum of 

AI { R \ 

— J 2 — - 2 10- 3 O 2 on a polar orbit of radius r ~ 7000 km. 
r \ r I 
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